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Abstract. The Ehrhart polynomial of a convex lattice polytope counts integer points in 
integral dilates of the polytope. We present new linear inequalities satisfied by the coeffi- 
cients of Ehrhart polynomials and relate them to known inequalities. We also investigate 
the roots of Ehrhart polynomials. We prove that for fixed d, there exists a bounded region 
of C containing all roots of Ehrhart polynomials of d-polytopes, and that all real roots of 
these polynomials lie in [— d, [d/2j). In contrast, we prove that when the dimension d is 
not fixed the positive real roots can be arbitrarily large. We finish with an experimental 
investigation of the Ehrhart polynomials of cyclic polytopcs and 0/1-polytopes. 



1. Introduction 

In this article, a lattice polytope P C M. d is a convex polytope whose vertices have integral 
coordinates. (For all notions regarding convex polytopes we refer to I n 1967 Eugene 

Ehrhart proved that the function which counts the lattice points in the n-fold dilated copy 
of P, 

i P : N -> N, i P (n) = # (nP n 7L d ) , 
is a polynomial in n (see [HIE] and the description in [Sj). In particular, ip can be naturally 
extended to all complex numbers n. In this paper we investigate linear inequalities satis- 
fied by the coefficients of Ehrhart polynomials and the distribution of the roots of Ehrhart 
polynomials in the complex plane. 

The coefficients of Ehrhart polynomials are very special. For example, it is well known that 
the leading term of ip(n) equals the volume of P, normalized with respect to the sublattice 
7L d H aff(P). The second term of ip(t) equals half the surface area of P normalized with 
respect to the sublattice on each facet of P, and the constant term equals 1. Moreover, the 
function i P (n) counting the number of interior lattice points in nP satisfies the reciprocity 
law ip(-n) = (-l) dimP i° P (n) HH HB| ■ 

Our first contribution is to establish new linear relations satisfied by the coefficients of 
all Ehrhart polynomials. This is a continuation of the pioneering work of Stanley, Betke 
& McMullen, and Hibi EH HI wno established several families of linear inequali- 
ties for the coefficients (see Theorems 13.11 and 13. 4j) . If we think of an Ehrhart polynomial 
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ip(n) = CdX d + Q_ix d_1 + • — h C\X + 1 as a point in d-space, given by the coefficient vector 
(q, Cd-i, . . . , Ci), their results imply that the Ehrhart polynomials of all <i-polytopes lie in a 
certain polyhedral complex. Betke and McMullen raised the issue page 262] of whether 
other linear inequalities are possible. We were indeed able to find such new inequalities in the 
form of bounds for the k-th difference of the Ehrhart polynomial ip(n). These are defined 
recursively via 

Aip(n) = ip{n + 1) — ip(n) 

and 

A k i P (n) = A (A k -H P (n)) for k > 1 and A°i P (n) = i P (n). 
Our first result (proved in Section^ is as follows. 

Theorem 1.1. If the lattice d-polytope Pci d has Ehrhart polynomial ip(n) = Cdn d + - ■ • + 
Co, then 

(^jA k i P (0) < ^JA^p(O) forO<k<£<d. 



In particular (put k = resp. £ = d), 

1) £ ^'<°< - u 



< A fc 2p(0) < ( d .)d\c d for 0<k<d. 



In Section |21 we give a proof of Theorem 11.11 using the language of rational generating func- 
tions as established in anc ^ ma ke a summary of known linear constraints and their 
strength. 

The relation between the coefficients and the roots of polynomials, via elementary symmetric 
functions, suggests that once we understand the size of the coefficients of Ehrhart polynomials 
we could predict the distribution of their roots in the complex plane. The second contribution 
of this paper is a general study of the roots of Ehrhart polynomials. 

There is clearly something special about the roots of Ehrhart polynomials. Take for instance 
the integer roots: Since a lattice polytope always contains some integer points (namely, its 
vertices), all integer roots of its Ehrhart polynomial are negative. More precisely, by the 
reciprocity law, the integer roots of an Ehrhart polynomial are those — n for which the open 
polytope nP° contains no lattice point. For instance, the Ehrhart polynomial { n ~^ d ) of the 
standard simplex in M. d (with vertices at the origin and the unit vectors on the coordinate 
axes) has integer roots at n — —d, —d + 1, . . . , — 1. 

The roots of the Ehrhart polynomial of the cross polytope 

O d = {(x 1 ,...,x d ) eR d : N + --- + |z d | < 1}, 
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also exhibit special behavior: Bump et al. [3] and Rodriguez [T§] proved that the zeros of i d 
all have real parts equal to —1/2. 

Using classical results from complex analysis and the linear inequalities of Theorem 13.51 we 
derive in Section 0] the following theorems: 

Theorem 1.2. (a) The roots of Ehrhart polynomials of lattice d-polytopes are bounded 
in norm by 1 + (d+ 1)! . 
(b) All real roots of Ehrhart polynomials of d-dimensional lattice polytopes lie in the 
half-open interval [—d, \_d/2\). 

The upper bound we present in Theorem II. 21 (b) is not tight. For example, in Proposition ^. 7[ 
we give a very short self-contained proof of the fact that Ehrhart polynomials for polytopes 
of dimension d < 4 have real roots in the interval [0, 1). In contrast with the above theorem 
we can also prove the following result. 

Theorem 1.3. For any positive real number t there exist an Ehrhart polynomial of suf- 
ficiently large degree with a real root strictly larger than t. In fact, for every d there 
is a d-dimensional 0/1-polytope whose Ehrhart polynomial has a real zero ad such that 
lim^oo a d /d = l/(2vre) = 0.0585 

Our third contribution is an experimental study of the roots and coefficients of Ehrhart 
polynomials of concrete families of lattice polytopes. Our investigations and conjectures 
are supported by computer experimentation using LattE |U E] and polymake jT2]. For the 
complex roots, we offer the following conjecture, based on experimental data. 

Conjecture 1.4. All roots a of Ehrhart polynomials of lattice d-polytopes satisfy —d < 
Re a < d — 1 . 

We also computed the Ehrhart polynomials of all 0/1-polytopes of dimension less than or 
equal to 4 and for many cyclic polytopes: 

Conjecture 1.5. For the cyclic polytope C(n, d) realized with integral vertices on the moment 
curve u d (t) := (t, t 2 , . . . , t d ) , 

ic(n,d)( m ) = vol(C(n,d))m d + i C ( n ,d~i)( m )- 

Equivalently, 

d 

ic{n,d){m) = >^ volfc (C(n, k)) m k . 
We have experimentally verified this conjecture in many cases. 
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2. AN APPETIZER: DIMENSION TWO 



Since Ehrhart polynomials of lattice 1-polytopes (segments) are of the form in + 1, where I 
is the length of the segment, we know everything about their coefficients and roots: the set 
of roots of these polynomials is {—l/£ : £ > 1} C [— 1, 0). 

The first interesting case is dimension d = 2. Pick's Theorem tells us that the Ehrhart 
polynomial of a lattice 2-polytope P is 

ip(n) = c 2 n 2 + cin + 1 , 

where c 2 is the area of P and c\ equals 1/2 times the number of boundary integer points of 
P. In 1976, Scott established the following linear relations. Two polytopes are unimodularly 
equivalent if there is a function which maps one to the other and which preserves the integer 
lattice. 

Theorem 2.1. ^7] Let ip(n) = c 2 n 2 + C\n + 1 be the Ehrhart 'polynomial of the lattice 
2-polytope P. If P contains an interior integer point, and P is not unimodularly equivalent 
to conv {(0, 0), (3, 0), (0, 3)}, then 

Cl < X -c 2 + 2 . 



By Pick's Theorem, for 2-polytopes with no interior lattice points, we have c\ — c 2 + 1. For 
P = conv {(0,0), (3,0), (0,3)}, we obtain i P (n) = 9/2 n 2 + 9/2n + 1. 

It is interesting to ask which degree-2 polynomials can possibly be Ehrhart polynomials. 
Since the constant term has to be 1, we can think of such a polynomial as a point (c 2 , c\) in 
the plane. From the geometry of lattice 2-polytopes, we know such an Ehrhart polynomial 
must have half-integral coordinates. Aside from Scott's inequality, we can trivially bound 
Cx > 3/2, since every lattice 2-polytope has at least 3 integral points, namely its vertices. 
From these considerations, we arrive at Figure ^ which shows regions of possible Ehrhart 
polynomials of 2-polytopes. 

Depicted are (parts of) three lines: 

(i) C! = 3/2 

(ii) ci = c 2 /2 + 2 

(iii) ci = c 2 + 1 
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Figure 1. Regions in which Ehrhart polynomials of lattice 2-polytopes lie. 
It consists of 3 half lines, an open region (only points with half-integral coor- 
dinates are possible), plus an exceptional point. 



and the point (c2,ci) = (9/2,9/2). The ray (0) shows the lower bound c\ > 3/2. This is 
a sharp lower bound, in the sense that we can have polygons with exactly three bound- 
ary integer points but arbitrarily large area. The ray (jnl) is Scott's bound, and the point 
(c2,Ci) = (9/2,9/2) corresponds to the "exceptional" polytope conv {(0, 0), (3, 0), (0, 3)} in 
Theorem 12.11 The rectangles conv {(0, 0), (2, 0), (2, x), (0, x)}, where x is a positive inte- 
ger, show that there is a point on (jnj) for every half integer. Finally, (jmj) corresponds to 
2-polytopes which contain no interior lattice point. There is a point on (jmj) for every half 
integer, corresponding to the triangles conv {(0, 0), (1, 0), (0, x)} for a positive integer x. The 
rays (0 and (|m|) meet in the point (1/2,3/2), which corresponds to the standard triangle 
conv {(0, 0), (1, 0), (0, 1)}. So the polyhedral complex containing all Ehrhart vectors consists 
of the polyhedron bounded by (0, (juj), and (00) (shaded in Figure |TJ), plus the ray (00), plus 
the extra point (02, C\) = (9/2, 9/2). In fact, only points with half-integral coordinates inside 
the complex are valid Ehrhart vectors. From these constraints, we can locate possible roots 
of Ehrhart polynomials of lattice 2-polytope fairly precisely. 
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Theorem 2.2. The roots of the Ehrhart polynomial of any lattice 2-polytope are contained 
in 



Proof. We consider three cases, according to Scott's Theorem 12.11 First, if the lattice 2- 
polytope P contains no interior lattice point then ip(n) = An 2 + (A + l)n + 1 (by Pick's 
Theorem), where A denotes the area of P. The roots of ip are at —1 and —1/A. Note that 
A is half integral. 

The second case is the "exceptional" polytope P = conv {(0, 0), (3, 0), (0, 3)} whose Ehrhart 
polynomial ip(n) = 9/2n 2 + 9/2n + 1 has roots —2/3 and —1/2. 

This leaves, as the last case, 2-polytopes which contain an interior lattice point and which 
are not unimodularly equivalent to conv {(0, 0), (3, 0), (0, 3)}. The corresponding Ehrhart 
polynomials ip(n) = c 2 n 2 + c\n + 1 satisfy the Scott inequality C\ < c 2 /2 + 2. Note that 
(because P has an interior lattice point) the area of P satisfies c 2 > 3/2. We have two 
possibilities: 

(A) The discriminant c\ — 4c 2 is negative. Then the real part of a root of ip equals — 
(which is negative). By Pick's Theorem c\ = c 2 — / + 1 where I is the number of interior 
lattice points, that is, — = — ^ — For fixed area c 2 , this fraction is minimized when 
I is smallest possible, that is I — 1. The imaginary part of a root of ip is plus or minus 



here we used ci > 3/2. As a function in c 2 , this upper bound is decreasing for c 2 > 1. Since 
c 2 > 3/2, we obtain as an upper bound for the magnitude of the imaginary part of a root 






(B) The discriminant c\ — 4c 2 is nonnegative. Then the smaller root of ip is 



— — — \Tc\--Ao2 > 

2c 2 2c 2 V 1 
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(Note that in this case c<i > 4.) 

Finally, the larger root is negative, since all the coefficients of ip are positive. □ 

3. Linear inequalities for the coefficients of Ehrhart polynomials 

In this section, we prove Theorem 1 which bounds the ratio of the fc-th and £-th differences 
of any Ehrhart polynomial solely in terms of d, k, and I. It is perhaps worth observing that 
most of our arguments are valid for a somewhat larger class of polynomials. To describe this 
class, we define the generating junction of the polynomial p as 

Sp(x) = J>(n) z B . 

n>0 

It is well known (see, e.g., [221 Chapter 4]) that, if p is of degree d, then S p is a rational 
function of the form 



(1) S p (x) 



[l - XY+ 1 



where / is a polynomial of degree at most d. Most of our results hold for polynomials p 
for which the numerator of S p has only nonnegative coefficients. Ehrhart polynomials are a 
particular case, as seen from the following theorem of Stanley. 

Theorem 3.1. |T5l Theorem 2.1] Suppose P is a convex lattice polytope. Then the generating 
function ^2 n>0 ip{n) x n can be written in the form of (fT| , where f(x) is a polynomial of degree 
at most d with nonnegative integer coefficients. 

Another well-known (and easy-to-prove) fact about rational generating functions (see, e.g., 
Chapter 4]) is the following. 



Lemma 3.2. Suppose that p G M[n] is a polynomial of degree d with generating function 
S p (x) = (adX d + • • • + a±x + «o)/(l — x) d+1 . Then p can be recovered 



(IS 



/ d + n — j 



d 



(2) p{n) = Y, a i 

3=0 

More generally, we have the identity 

(3) A k p(n) = J2aJ d + "7 j ) fork>0. 



3=0 
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Proof. Equation (J2J) follows from expanding 1/(1 — x) d+1 into a binomial series. For (|3*|). we 
proceed by induction on k. For = 0, the statement is (0), while for k > 1 we have by the 
induction hypothesis 

A k p(n) = A fc_1 p(n + 1) - A fc_1 p(n) 

E//d + n+ l — A fd + n — j 
aj 

3=0 
d 

= ^2 a i 

3=0 

□ 



d - A; + 1 J \d - k+1 



d + n — j 
d — k 



Combining Theorem 13. II and Lemma fa .21 immediately yields the following fact. 
Corollary 3.3. For any lattice polytope P and k>0, we have A fc zp(0) > 0. 

Proof. This follows because those binomial coefficients in the final expression for A h p(n) are 
either positive or zero. □ 

Proof of Theorem \l.l\ We will use the falling-power notation dt = d(d — 1) ■ • • (d — j + 1), 
along with the obvious relation k?~ < B- for j < k < £, and the identity 

'd - j\ fd\ ki 
K d-k) \k) dl' 
The statement now follows from Lemma 13.21 (j3J) by 

d\ fd - A /d\ (d\ kl /d\ /d\ £i fd\ (d - j 

i)\d-k) ~ \e)\k)dl < \k)yi)a% ~ \k)\d-£, 

□ 

Theorem 11.11 is not the first set of linear inequalities on coefficient vectors of Ehrhart poly- 
nomials. Indeed, in 1984, Betke and McMullen P3 Theorem 6] obtained the following in- 
equalities. 



Theorem 3.4. Let P be a lattice d-polytope whose Ehrhart polynomial is Y2t=o c i nl - Th 



CI) 



c r < (-l) d - r s(d, r) c d + (-1)^-1 <*> r + [) for r = l,2,...,d-l, 
where s(k,j) denote the Stirling numbers of the first kind. □ 
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In that paper, Betke and McMullen sent out a challenge to the community to discover new 
inequalities for these coefficient vectors. The following theorem sums up the current state of 
affairs. 



Theorem 3.5. Let P be a d- dimensional lattice polytope, with Ehrhart polynomial ip{n) = 

Y2i=o c i n% = Y2i=o a i( U+ d~ l ) ' Then ^ e following inequalities are valid for < k < £ < d and 
0<i<d: 



(4) 



C r < 



\d—r 



s(d,r) c d + 



\d-r-l 



s{d,r+l) 

(d-iy- 



(5) (T)AV(O) > (f)A fc z P (0), 



d + 1 



Cd > c d - 



(6) 
(7) 

(8) A k i P (0) > 



(9) 
(10) 



c d > c /d\, 

d+1 

Cd-l > c 



2(d-l) 



ip(l) > d+1, 
k. 



'12) cu > 0. 



Moreover, 

(13) + arf_i + ■ ■ ■ + a d _i < a + a\ + ■ ■ ■ + a,i + a i+1 for all < i < [(d — 1)/2J 



Whenever a s ^ £m£ a s+ i = • • • = a d — 0, then 

(14) a + ai H hft; < a s + a s _i H h a s _j 

finally, if a d ^ 0, i/ien 

(15) ai < a. 



/or a// < z < s ; 
/or a// 2 < % < d . 



Proof. The inequalities (jlj) and (JHJ) are the contents of Theorems 13.41 and II. 11 while (JBJ), (JZJ), 
and (JHJ) are the special cases (A;, £) = (d — (A;, £) = (0,1), and k — 0, respectively. 

(J2J) and (fTUj) say that the volume and the normalized surface are at least as big as for a 
primitive simplex. Inequality (fTTj) follows from Ehrhart reciprocity. Inequality ()12j) is the 
statement of Theorem 13 .11 Incidentally, Q also follows from (JSJ), and (fTTj) follows from ()12J) . 
both by specializing to i = d. Inequality (fT4"j) was proved by Stanley [21], and inequalities 
(fT3J) and (HSJ) by Hibi HUGH]. □ 



It is illuminating to compare these inequalities with each other. Since inequality (|12j) was 
used to prove Theorem 13.41 (by Betke and McMullen) and Theorem 11.11 it seems stronger 
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than the other inequalities. Indeed, the only inequality among (|4"])-([12j) which does not follow 
from (jl2j) is ()10|). Experimental data for small d shows that neither Q nor (jSJ) imply the 
other. 

The set of linear inequalities of Theorem 13.51 describes an unbounded complex of half-open 
polyhedra in inside which all coefficient vectors of Ehrhart polynomials live. From 

this, we obtain a bounded complex Q d by cutting with the normalizing hyperplane q = 1. 
By (|14J1 each constraint a s 7^ 0, a s+ i = • • • = = for s = 1,2, ... ,d defines a half-open 
polytope E s G Q d of dimension s that is missing one facet; E is a single point. 

Here are some particular cases: The bounded complex Q 3 consists of one half-open s- 
dimensional simplex for each s — 0, 1, 2, 3 (Figure |2J), and the half-open 3- and 4-dimensional 
polytopes of Q 4 are shown in Figure El 




Figure 2. The complex Q 3 of half-open polytopes, inside which the possible 
Ehrhart coefficients of all 3-dimensional polytopes lie. The facets of the tetra- 
hedron corresponding to a 3 ^ are a > 0, c 2 > 1, (|T4*jl and (fT3|) ; those of the 
triangle corresponding to 03 = 0, a 2 7^ are ao, «i > and (JHJ); ari d those of 
the segment a 2 = 03 = 0, a\ 7^ are ao > and (JHJ). 

An important question about any linear inequality is whether or not it defines a facet of Q d . 
We rephrase Betke and McMullen's question 0: 
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Figure 3. The 4-dimensional (top) and 3-dimensional (bottom) member of 
the complex Q 4 . The facets of the 4-dimensional polytope are a > 0, C3 > 
5/12, (H1J) for % = 0, (JT3J for % = 0,1, and JTSJ) for i = 2,3; those of the 3- 
dimensional one corresponding to 04 = but 03 7^ are a ,ai > 0, C3 > 5/12, 
(O for i = 1, (HU) for i = 1, and (fT5j) for i = 2; etc. 

Problem. Are there other linear inequalities for the coefficients of an Ehrhart polynomial 
aside from those in Theorem \3.bl ? Do they define facets of the polyhedral complex inside 
which all coefficient vectors of Ehrhart polynomials live? 



4. The roots of Ehrhart polynomials 

When one has a family of polynomials, a natural thing to look at are its roots. What 
is the general behavior of complex roots of Ehrhart polynomials? As a consequence of the 
inequalities on its coefficients, we give bounds on the norm of roots of any Ehrhart polynomial 
in dimension d. The basis { ~ J ) : < j < d} of the vector space of polynomials of degree d 
turns out to be much more natural than the basis {n l : < i < d} for deriving bounds on the 
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roots of Ehrhart polynomials ip(n) = Yli=o ai { n+ d~ l ) = Sf=o c * nl - Also, recall the following 
classical result of Cauchy (see, for example, [E3 Chapter VII]). 

Lemma 4.1. The roots of the polynomial p(n) = Cdn d + Cd-\n 1 + • • • + Co lie in the open 
disc 

{z G C : \z\ < 1 + max 
o<i<d 

□ 




Now we study roots of Ehrhart polynomials in general dimension. We first give an easy proof 
bounding the norm of all roots. 



Proof of Theorem M.^f a). By Lemma 14.11 and Theorem 13.41 the maximal norm of the roots 
of ip is bounded by 



1 + max 

0<j<d 



Cd 



< 1 + max 

0<j<d 



(-iy->s(d,j) + (-ir- ^\ +1) 

[a - l)\c d 



< l + d\ + d\d = l + (d+l)\. 
Here we have used the estimate s(d,j) < \s(d,j)\ < d\ and the fact that Cd > l/d\. 



□ 



While using crude estimates gives us a bound of 1 + (d+ 1)!, which makes the main point that 
there exists a bound dependent only on d, the actual bound on the roots can be improved 
greatly for specific values of d. First of all, for small d, we can compute the inequalities 
exactly; here the inequalities from Theorem 11.11 are used along with the Betke-McMullen 
inequalities. This gives appropriate bounds on the ratios of the coefficients of the Ehrhart 
polynomial. Second of all, Lemma 14.11 is not the best tool to use for specific cases, since 
calculating the inequalities for small d yields much lower bounds for Ci/c d when i is large. 
Instead, we use the following proposition. 

Proposition 4.2 (Theorem 27.1 ^Sj). Let p(n) = Cdn d + Cd-\n d ~ 1 + ■ ■ - + co be a polynomial. 
Then the maximal value of the norm of a root of p(n) is the value of the maximal root of 
p'(n) = \c d \n d - {cd-iln 11 ' 1 - \c d -2\n d ~ 2 \c \. □ 



We use this and the exact calculation of the inequalities in question to obtain the following 
tighter bounds on the roots of Ehrhart polynomials of d-polytopes. 



d 


2 


3 


4 


5 


6 


7 


8 


9 


bound 


3.6 


8.5 


15.8 


25.7 


38.3 


53.5 


71.4 


92.0 
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The bound appears to grow roughly quadratically. We suspect that there is a bound for the 
roots of Ehrhart polynomials of G?-polytopes which is polynomial in d. For real roots this is 
certainly the case; we prove next that all real roots of Ehrhart polynomials of d-polytopes 
lie in the interval [—d, \_d/2\). For this, we will use the following well-known bound. 

Lemma 4.3. (Newton Bound) Let f e K[n] be a polynomial of degree d and B e K be such 
that all derivatives of f are positive at B : f^\B) > for £ — 0, 1, . . . , d. Then all real roots 
of f are contained in (— oo,B). □ 



Proof of Theorem \l.^f b). The lower bound follows from Theorem 13 . 1 1 and the simple obser- 
vation that for (real numbers) n < —d the binomial coefficients in 

are all positive or all negative, depending on the parity of d. 

As for the upper bound, let B = [_d/2\. We now show that a < B for any real root a 
of ip(n). For this, we will make use of the fact that the second highest coefficient of any 
Ehrhart polynomial measures half the normalized surface area. This coefficient reads 

1 d 

Cd ~ l = (d-iv ^ ai ^~ 2z + 1 ) 

^ i=0 

when expressed in terms of the Oj's, so that the following inequality is valid: 

d 

(16) (d-iy.ca-! = ^ ai (d-2i + l) > 0. 

Note that the coefficient s(i) = d — 2i + 1 of a\ in (|16j) is positive for < % < B and non- 
positive for B + 1 < % < d. We now express the £-th derivative of ip evaluated at n = B as 
ip\B) = (£\/d\) J^f =0 a>i gi(B, £), and claim that for < £ < d, there exists a X(£) > with 

gi (B, t) > X(£) s(i) for all < i < d. 

This claim is the statement of Lemma 14.51 below. The proof of Theorem ll.2f b) now follows 
from this relation, inequality (J16)) . ao = 1 and > for 1 < i < d via the following chain of 
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inequalities: 



d 



< 



i=0 



d 



d 



< 



\(£)s(i)) ai + \(£) ^s(i)oi 



i=0 



i=0 



d 



J29i(B,£) a, 




i=0 



□ 



Remark. It is a well-known fact that Ehrhart polynomials of lattice polytopes form a special 
class of Hilbert polynomials. More strongly, they are special examples of Hilbert polynomials 
of Cohen-Macaulay semi-standard graded fc-algebras [25 (this is essentially the content of 
Theorem IH.ljl . It is then natural to ask whether Ehrhart polynomials are special or whether 
the bounds proved above hold in more generality We stress that inequality (jlfij). used in 
previous arguments, comes from geometric information about Ehrhart polynomials ip(n). 
Indeed, from the following proposition and Theorem II. 2f b). Ehrhart polynomials are special 
in their root distribution: 

Proposition 4.4. For degree d Hilbert polynomials associated to arbitrary semi-standard 
graded k-algebras the negative real roots are arbitrarily small and d—1 may appear as a root. 
In contrast, for fixed degree d, Hilbert polynomials of Cohen-Macaulay semi-standard graded 
k-algebras have all its real roots in the interval [—d, d — 1). 

Proof. Indeed, it follows from Theorem 3.8] that for fixed d and positive integers a , . . . , a^, 
the polynomial Oq(x + a\)(x + a 2 ) . . . (x + a a) is the Hilbert polynomial of a semi-standard 
graded /c-algebra. Also, observe that the chromatic polynomial of the complete graph on 
d vertices has highest root d—1, and that chromatic polynomials are known to be Hilbert 
polynomials of standard graded algebras by a result attributed to Almkvist (see the proof 
given by Steingrimsson j2H])- Thus the first statement holds. 

Now, in a Cohen-Macaulay semi-standard graded algebra, the Hilbert polynomial can be 
written as p(n) = Ylt=o a «( n+ d _ *)' where a« > for < i < d. Observe that all the binomial 
coefficients in p(n) are positive for (real numbers) n > d — 1, which establishes the upper 
bound of d — 1. For the lower bound, observe that for (real numbers) n < —d all the binomial 
coefficients are positive, respectively negative, depending on the parity of d. □ 

To complete the proof of Theorem 1.2 (b), we need only to prove the following lemma. 
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Lemma 4.5. Fix < £ < d — 1 and consider again the functions s,g : {0,1, ... ,d} — > Z 
defined by s(i) = d — 2i + 1 and g{i) = gi(B, £). Moreover, if we set 

(17) X(i) = l(g(B)-g(B + l)) = d - £ 11^-*) > °> 

rp ( [d-i] \ kei 

then 

(18) g(i) > \(£)s(i) for 2 = 0,1,..., d. 

For this, we will also need to prove Lemma l4~B1 b elow . We will write [d— l]o = {0, 1, . . . , d— 1}, 
[d— 1] = {1, 2, . . . , d— 1}, and ( ) for the set of all t-element subsets of the finite set S. Now 
we express ip(n) as 

d d-l 



ip(n) = ^ Qj 1~[ (n + d - i - k), 

i=0 k=0 

so that the £-th derivative of ip is 

£\ d 

= dJ.^ ai II (n + d-i-k) 



o Ie ad-i] \ ke[d-l] \i 
J\12 ai H ]J(n + d-i-k) 



Note that we now have an explicit formula for the coefficient of a, in (d\/£\) ip. 

(19) 9l (n,£) = \{{n + d-i-k). 

The following lemma shows that the piece- wise linear function interpolating g : {0, 1, . . . , d} — > 
Z, g{%) = gi(B,£) is positive, and its slope weakly increases in the range < £ < d and 
< i < B + 1. See FigureH 

Lemma 4.6. The following inequalities are satisfied for < £ < d and < % < B + 1: 

(20) 9i (B,£) > 0, 

(21) 9t{B, £) — gi+i(B, £) > g i+1 {B, £) - g t+2 (B, £). 
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Proof. Equation (|20|) follows because k < d — 1 and % < B + 1 imply B + d — i — k > 0. To 
show (J21J) . we abbreviate m := B + d — i and inspect the difference 

gi(B,e)-g i+1 (B,t) = Yl \[(rn-k)- ^ J] (m - fc - 1) . 

H [d i-l°) ka M ld d --l°) keJ 

If ^ J, then the term corresponding to / in the first sum cancels with the term corresponding 
to J = {i — 1 : i G /} in the second sum: 

Y[(m-k)-Y[{m-k-l) = JJ ((m - k) - (m - (k - 1) - l)) = 0, 

kei fceJ kei 

so we are left with summing over the sets / G (^T 1 ] ) that contain and the sets J that 
contain d — 1. But for such summation sets, the difference simplifies to 

gi (B,£)-g i+1 (B,£) = (m-0) ]T J] (m - A;) - (m - d) ^ J] (m - (k + 1)) 



d J(B-d-z-ife) 



^| fee/ 

and (j2*Tj) follows by comparing the expressions gi(B, £) —g i+ i(B, £) and g i+ i(B, £) —gi +2 (B, £) 
term by term. □ 

In the following, we will use Iverson's notation (see the expression [S] evaluates to 1 
resp. according to the truth or falsity of the logical statement S. 

Proof of Lemma \4 ■ 5\ First note that s(B) = 1 for even d, so that 

g(B)-Xs(B) = l g (B) + lg(B + l) > 0; 

for odd d, we have s(B + 1) = 0. Now note that the graph of (the piecewise-linear function 
interpolating) As is a line, while g(B + 1) > by (|2(Jj) and the slope of the graph of g is 
weakly increasing on [0, B + 2] by (|2T|) (see FigureEJ); this proves (JTHJ) for < i < B+ [d odd]. 

Set j = i — B, so that we still need to prove (fTHj) for 1 + [d odd] < j < d — B. By plugging 
(1191) and ()17|) into (|18|) and rearranging, we must show that for these values of j 

( 22 ) e n^-*-*) + ~ (2.7 - odd] - 1) £ n^-^) > °- 



Note that each term in the second sum of ()22|) is positive, and decompose the index sets / 
in the first sum into disjoint unions I = I + U K such that I + C {0, 1, . . . , d — 2j} and 
K C {d — 2j + 1, . . . , d — 1}, and therefore d — j — k > for all k G i+. 
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Figure 4. The graphs of the functions g and As (solid for odd d, dashed for even d). 
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If |7T| is odd, then the summand cr(Jl") corresponding to I + U K cancels with the one corre- 
sponding to J + U (d — j — K), so we only need to consider even \K\. In that case, o~(K) > 
(resp. cr(K) < 0) if \K D [d — j + 1, d — 1]| is even (resp. odd). In total, there are more than 
enough positive terms in f)22|) to cancel the negative summands. □ 

Proposition 4.7. VFe /mve a < 1 for any real root a of an Ehrhart polynomial ip of a 
lattice polytope P of dimension d < 4 . 

Proof. It is enough to prove the statement in dimension 4 because of Theorem II. 2f b). Sup- 
pose f(n) = pn A + qn 3 + rn 2 + sn + 1 is the Ehrhart polynomial of a lattice 4-polytope P. 
We know p > and q > 0. Because /(l) counts the lattice points in P, we know that 
p + q + r + s + 1 > 5. By the reciprocity law, /(— 1) > 0, so p — q + r — s + 1 > 0. The top two 
coefficients of the shifted polynomial g(n) — f(n + 1) = pn 4 + (4p + g)n 3 + g2n 2 + + go 
are positive, as is the constant term g = g(0) = /(l). We will show that g 2 and g\ are 
nonnegative, and hence, by Descartes' rule of signs, g does not have a positive root. This 
implies that fin) — gin — 1) does not have a real root larger than 1. To prove that gi > 0, 
we add the inequalities /(l) > 5 and /(— 1) > to obtain 2p + 2r > 3 or r > | — p, whence 
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g 2 = 6p + 3g + r > 5p + 3q + | > (because p,q > 0). A similar reasoning yields 

3 

gi = 4 p + Sq + 2r + s = {p + q + r + s) + (3p + 2g + r) > 4 + 2p + 2g + - > ; 
here we used the inequality /(l) > 5 again. □ 

We now conclude with the proof of Theorem 11.31 

Proof of Theorem M.iA Given a positive integer d, consider the convex polytope Pd defined 
by the facet inequalities: 

< x < x k < 1 for 1 < k < d - 1. 

Pd is an order polytope in the sense of [20] and thus it has 0/1 vertices. We claim that the 
Ehrhart polynomial of Pd is given by ip d (n) = (Bd(n + 2) — Bd(0))/d where Bd(x) is the d-th 
Bernoulli polynomial. 

Indeed, from the facet-defining inequalities of Pd one sees that ip d (n) is the number of d- 
tuples of nonnegative integers (a , a±, . . . , a^-i) such that a < a k < n. If a = j then 
there are n — j + 1 choices for each a k (k > 0). Hence ip d (n) = X]j=o( ri — + l) <i_1 = 
Si=i J d l - classical identity of Bernoulli that says X]fc=o = (Bd(n) — Bd(0))/d. Thus 
we get ip d {n) = (B d (n + 2) — B d (0))/d, a polynomial of degree d. Note that when d is 
odd then 5^(0) = 0. Finally, the results of [21] imply that the largest real zero of Bd(n) is 
asymptotically d/(27re). Therefore, as stated, as the degree d grows, the Ehrhart polynomial 
of Pd has larger and larger real roots. It is worth remarking that since d is the degree of the 
Ehrhart polynomial of Pd it differs from the upper bound of \_d/2\ in Theorem 11.21 only by 
a constant factor. □ 



5. Special families of polytopes 



We begin this section with some charts showing the behavior of roots for hundreds of Ehrhart 
polynomials computed using LattE and Polymake. In Figure El we show the distribution of 
roots of a large sample of Ehrhart polynomials of lattice 3-polytopes. 

In Table we collected a small sample of Ehrhart polynomials of 0/1 polytopes and cyclic 
polytopes from the experiments we performed. 



Ehrhart polynomial P(s) 



cube 



t A + 3t 2 + 3t + 1 



cube minus corner 



5/6t 3 + 5/2t 2 + 8/3* + 1 
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prism 



nameless 



2/3t 3 + 2t 2 +7/34 + 1 



2/3 + 2i 2 + 7/3t + 1 



octahedron 
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bypyramid 



1/2 f 3 + 3/2t 2 + 2t + 1 



1/6 f 3 + f 2 + gj_+ 1 



unimodular tetrahedron 



fat tetrahedron 
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cut(4) 



■4/3x 3 



cyclicOl: 5-8. poly 
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15 



. x 4 + 5/4x 3 + 



~7TT 



1 



halfcube(5) 



+ 11/3 x 4 + 16/3 x 3 + 10/3 x 2 + 9/5 x + 1 



Cyclic(2,5) 



10x 2 



- 4x 



1 



Cyclic(3,5) 



16x 3 



10x 2 



- 4x + 1 



Cyclic(4,5) 



12x 4 



16x 3 



10x 2 



- 4x + 1 



Cyclic(2,6) 



20x 2 + 5x + 1 



Cyclic(3,6) 



70x 3 + 20x 2 + 5x + 1 



192x 4 + 70x 3 + 20x 2 + 5x + 1 



Cyclic(4,6) 



288x 5 + 192x 4 + 70x 3 + 20x 2 + 5x + 1 



Cyclic(5,6) 



Cyclic(2,7) 



35x 2 + 6x + 1 



224x 3 + 35x 2 + 6x + 1 



Cyclic(3,7) 



1512x 4 + 224x 3 + 35x 2 + 6x + 1 



Cyclic(4,7) 



Cyclic(2,8) 



56x 2 + 7x + 1 



Cyclic(3,8) 



588x 3 + 56x 2 + 7x + 1 



8064x 4 + 588x 3 + 56x 2 + 7x + 1 



Cyclic(4,8) 



Table 1. The Ehrhart polynomials for some well-known lattice polytopes. The choice of 
coordinates for cyclic polytopes was t = 1, n. The rest are listed Ehrhart polynomials comes 
from 0/1 polytopes selected from Ziegler's list. It includes the Ehrhart polynomials of all 
3-dimensional 0/1-polytopes. 
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FIGURE 5. The zeros of Ehrhart polynomials corresponding to 100000 ran- 
dom 3-dimensional lattice simplices. 

5.1. 0/1-polytopes. We computed the Ehrhart polynomials for all 0/1 polytopes of dimen- 
sion less or equal to 4 (up to symmetry there are 354 different 4-polytopes) . In Figure El 
we plotted their roots. In our computations we relied on the on-line data sets of 0/1 poly- 
topes available from Polymake's web page and those discussed in Ziegler's lectures on 0/1 
polytopes [T3]. Several phenomena are evident from the data we collected. For example, 
in Table Q we see two combinatorially different polytopes that have the same Ehrhart poly- 
nomial. These are the so called "nameless" polytope of coordinates (1,0,0,0), (1,1,0,0), 
(1,0,1,0), (1,1,1,0), (1,0,1,1), (1,1,0,1) and the octahedron. Another example of regu- 
lar distribution appears also in Figure |3 We show the roots of the Ehrhart polynomials 
associated to the Birkhoff polytope of doubly stochastic n x n matrices for n = 2, . . . , 9. 

5.2. Cyclic polytopes. Cyclic polytopes form a family whose combinatorial structure (i.e. 
/-vector, face lattice, etc) is well understood. The canonical choice of coordinates is given 
using the moment curve 



A cyclic polytope is obtained as the convex hull of n points along the moment curve. Thus 
we fix ti, t 2 , ■ ■ ■ , t n and define C(n, d) := conv{v d (ti), v^fo), ■ ■ ■ , Vditn)}- Cyclic polytopes 



(23) 
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are lattice polytopes exactly when ti G Z. There is a natural linear projection connecting 
these cyclic polytopes. 

Lemma 5.1. Consider the projection 7r : M. d — > that forgets the last coordinate. The 

inverse image under it of a lattice point y G C(n, d — 1) n Z^" 1 is a line that intersects the 
boundary ofC(n,d) in exactly two integral points. 



Proof. We need to prove that, given t±, t 2 , ■ ■ ■ , t d G Z and A 1; A 2 , . . • , A rf G 

d d 



fe=i 



fc=i 



For 1 < j < d, let yj = Ylt=i we know that y 1 , y 2 , . . . , G Z. We need to prove that 



y = {Vi, ■ ■ ■ , Vd) = ^2 Afe E ^ 



fc=i 



This identity means that y lies on the hyperplane spanned by Vd{ti), ■ ■ ■ , Vd{td), which can 
be expressed via a determinant: 



det 



1 1 

y Vd(ti) 



1 

Vd{td) 







Writing this determinant out through the first column and solving for yd gives 
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This expression yields an integer if we can prove that D divides the determinants appearing 
in the numerators. Equivalently, the substitution tj = tk in any of the numerators evaluates 
the determinant to zero, which is apparent. □ 

Consequently, Conjecture 11.51 is equivalent to saying that the number of lattice points in a 
dilation of a cyclic polytope by a positive integer m is equal to its volume plus the number 
of lattice points in its lower envelope. From the above lemma and Pick's theorem, it follows 
that Conjecture 11.51 is true for d = 2. 
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